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We extend the Coulomb and Proximity Potential Model (CPPM) to study two-proton (2p) radioactivity from 
excited states while the proximity potential is chosen as AW95 proposed by Aage Withner in 1995. Demon- 
stration reveals that the theoretical results acquired by CPPM exhibit a high level of consistency with prior 
theoretical models, such as the unified fission model (UFM), generalized liquid drop model (GLDM) and ef- 
fective liquid drop model (ELDM). Furthermore, within the CPPM, we predicted the half-lives of potential 
2p radioactive nuclei for which experimental data are currently unavailable. The predicted results were then 
assessed, compared with UFM, ELDM, and GLDM models, and examined in detail. 
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I. INTRODUCTION 


Except for the well-known a, (, y and cluster decays, 
etc[1—10], there exist some exotic modes of radioactivity in 
proton-rich nuclei [11—13], such as proton and two-proton 
(2p) radioactivity [14-18]. 2p radioactivity is an extremely 
exotic mode that is energetically possible in less massive nu- 
clei near the proton-drip line [19-25], and it can reveal abun- 
dant information of nuclear structure, including the nuclear 
radius, wave function of the emitted two protons, spin and 
parity, deformation effect and so on [26-28]. This novel de- 
cay mode was firstly predicted by Zel’dovich and Goldan- 
sky in 1960s [29-31]. Based on the pioneering work of 
Zel'dovich and Goldansky, an extensive range of theoretical 
models have been proposed to describe this exotic decay pro- 
cess [32-35]. However, due to the limitations in radioactive 
beam facilities and detection technology, 2p radioactivity was 
not experimentally confirmed until the observation of 4°Fe 
— Cr+p+p decay by Giovinazzo et al. at GANIL (France) 
and independently by Pfiitzner et al. at GSI (Germany) in 
2002, which provided the first experimental evidence of this 
decay mode [36, 37]. Since then 2p radioactivity has been 
recognized as the significant decay mode for proton-rich nu- 
clei [38], and it has been detected and studied in several nu- 
clei, such as the resonant ground state of ®Be [39, 40] and 
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120 [41], and in the excited state selectively populated in °Be 
[40] and !40 [42], etc. 


Numerous studies have shown that the 2p radioactivity not 
only occurs from ground state but also from short-lived ex- 
cited state. Janecke was the first to discuss the possibility of 
B-delayed 2p (82p) radioactivity [38], while Goldansky pre- 
dicted the occurrence of {2p radioactivity could be found in 
62p emitters of Z = 10 ~ 20. In 1983, Cable et al. re- 
ported the first experimental observation of 682p radioactivity 
[43]. Subsequently, an increasing number of (2p emitters 
have been detected using silicon detector telescopes, making 
it possible to measure the energy of two individual protons 
with high precision. Since Cable et al. [43] discovered the 
B2p radioactivity from ??Al, shortly they observed more 62p 
radioactivity from ?°P [44] and °°Ca [45]. Hereafter, several 
other 682p nuclei were found, including ?Si [46], 27S [47], 
31Ar [48], 3°Ti [49], 43Cr [50, 51] and °°Ni [52]. In addi- 
tion to populating excited state 2p radioactivity via 8 decays, 
2p radioactivity has also been observed from excited state fed 
by nuclear reactions such as pick up, transfer or fragmenta- 
tion, including 140 [42], '7!8Ne [53-57], Mg [58, 59] and 
28,295 [60, 61]. In 2006, Mukha et al. first reported 2p ra- 
dioactivity from °4Ag in an experiment at GSI [62]. 


From a theoretical perspective, several methods have been 
proposed over the past few decades to study the mechanism of 
2p radioactivity, including both microscopic and phenomeno- 
logical models over the previous decades. In general, there 
are three distinct ways for proton-rich nuclei to emit two pro- 
tons: (1) two-body sequential emission, (2) three-body si- 
multaneous emission, and (3) diproton emission (also called 
?He cluster emission). The *He cluster emission is an ex- 
treme scenario with the two strongly correlated protons that 
can only survive for a brief time before splitting after passing 
through the Coulomb barrier. Recently, based on the CPPM, 
Yao et al. [64] and Ghodsi et al. [65], as well as Deng et 
al. [66] and Santhosh et al. [67], performed comparative 


studies of various proximity potential formalisms to study a 
decay, proton radioactivity and cluster radioactivity, respec- 
tively. Considering that the 2p radioactivity decay process 
shares the same theory as a decay and proton radioactivity, 
i.e., barrier penetration [68-73], we extended the CPPM pre- 
sented in Ref. [74] to account for ground state 2p radioac- 
tivity. It was found that both the calculated and predicted 
results were highly consistent with experimental data and re- 
sults from other theoretical models. In fact, except for theo- 
retical models the empirical formulae are excellent tools for 
investigating 2p radioactivity involving a four-parameter em- 
pirical formula provided by Sreeja et al. [34] and a two- 
parameter empirical formula proposed by Liu ef al. [63]. 
Consequently, a desirable question is whether 2p radioactiv- 
ity from excited states can be regarded as *He cluster emis- 
sion and described by the CPPM. To address this question, we 
systematically studied the half-lives of 2p radioactivity from 
excited states within the CPPM for nuclei with 8 < Z < 47 
in this work. 

This paper is organized as follows: In Sec. II, the theo- 
retical framework of 2p radioactivity half-life and the CPPM 
formalism are described briefly. In Sec. III, detailed calcu- 
lations and discussion are presented. Finally, a summary is 
given in Sec. IV. 


Il. THEORETICAL FRAMEWORK 


85 
A. The 2p radioactivity half-life formalism 
The half-life of 2p radioactivity is generally determined by 
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89 (1) 
where À is the decay constant. v is the assault frequency re- 
lated to the harmonic oscillation frequency presented in the 
Nilsson potential [75], it can be written as 


41 
ae 
Here h, h, A and w are the Planck constant, reduced 
Plank constant, mass number of parent nucleus and angu- 
lar frequency, respectively. $2, G?[A/(A — 2))?"x? 
represents the preformation probability of the emitted two 
protons in parent nucleus obtained by the cluster overlap 
2n)! [20, 86] with n 


22" (nl)? 
(3Z Me 3 — 1 is the average principal proton oscillator quan- 
tum number and x? is set as 0.0143 acrooding to Ref. [23]. 
The penetration probability P can be calculated by the WKB 
approximation and written as 


hv = hw ~ (2) 
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where K(r) A) — Q2p|, pe denotes the reduced 


108 Mass [35]. r is the mass center distance between the daugh- 
107 ter nucleus and emitted two protons. fin and rout are classical 
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inner and outer turning points which can be obtained from the 
conditions V (rin) = V (Tout) = Q2p. The whole interaction 
potential V (r) between the emitted two protons and daughter 
nucleus is composed of the Coulomb potential Vo (r), nuclear 
potential Vy (r) and centrifugal potential Vp(7). It can be ex- 
pressed as 

V(r) = Vn (r) + Ve(r) + V(r). (4) 
In CPPM, the nuclear potential is replaced by proximity po- 
tential. A detailed description of the proximity potential will 
be provided in Sec. IIB. 

Assuming a homogeneous spherical charge distribution for 
the daughter nucleus, the Coulomb potential Vo(r) is postu- 
lated to be the potential of a uniformly charged sphere with 
radius R. It is expressed as follows: 


v= 


where R = Rı + Re is the separation radius with R, and Rə 
being the radii of daughter nucleus and emitted two protons, 
respectively. 

For V(r), we choose the Langer modified form, because 
L(+ 1) > (L+ $)? is a necessary correction for one- 
dimensional problem [76]. It can be written as 


ZZee" 13 (2), r < R, 


r>R, K 


ZıZ2e? 
r 7 


W(é+ 3)? 


V(r) = 2ur? 


; (6) 
where £ is the obrital angular momentum taken away by the 
emitted two protons. 


B. The proximity potential formalism 


The phenomenological proximity potential was first pro- 
posed by Blocki et al. [77] in 1970s for heavy-ion reactions 
in 1970s for the first time. It provides a simple formula for the 
nucleus-nucleus interaction energy as a function of the sepa- 
ration between the surfaces of the approaching nuclei, with 
adjustable parameters that makes use of the measured values 
of the nuclear surface tension and surface diffuseness. In the 
CPPM, the potential energy barrier is modeled as the sum of 
Coulomb potential, proximity potential and centrifugal po- 
tential for both the touching configuration and separated frag- 
ments. In this work, we select the AW95 proposed by Aage 
Withner in 1995 [78] as an example to replace nuclear po- 
tential. This set of proximity potential can be expressed as 


Vo 
V(r) = aa (7) 
1+ exp (Fes) 
Here Vo = 16r p4p- ya with 
Mp 1 fm. (8) 
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R; = (1.2A;/ — 0.09) fm (i = 1,2) Ro = Ry + Ro. The 
surface energy constant ~y has the form as 
)| » © 
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A 
where the coefficients yọ = 0.95 Mev/fm? and K, = 1.8. Here 
Ai, Zi and N; (i = 1, 2) are the mass number, proton number 
and neutron number daughter nucleus and emitted two pro- 
tons, respectively. 


Ill. RESULTS AND DISCUSSION 


Blocki first presented the proximity potential in 1977 to 
characterize the interaction potential between any two nuclei 
in the separation degree of freedom, based on the proximity 
force theorem [77]. Hence, numerous nuclear proximity po- 
tentials have been widely used in nuclear physics research 
[79-83]. The CPPM is a phenomenological model that is 
commonly used to study the two-body problem in the context 
of 2p radioactivity. It considers the interaction potential be- 
tween the parent and daughter nuclei as the sum of the nuclear 
potential, Coulomb potential, and centrifugal potential. Un- 
like other models, the CPPM utilizes the proximity potential, 
which replaces the nuclear potential with a simplified formal- 
ism based on the proximity force. This provides the model 
with the advantage of adjustable parameters, making it sim- 
ple yet accurate. In our previous work [74], we extended the 
CPPM to study 2p radioactivity from the ground state, the 
main intention of this work is to further extend the CPPM to 
investigate 2p radioactivity from excited states. 
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Fig. 1. (color online) The whole interaction potential of 77Mg* and 
29 qa* 
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Firstly, we performed the calculations on 2p radioactivity 
half-lives for 140*, '7Ne*, 18Ne*, ??Mg*, 29S* and 94Ag* 
(* represent the excited state), all of the calculated results are 
listed in the Table 1. In this table, the former four columns 
represent the 2p decay process, experimental 2p radioactiv- 
ity released energy, spin and parity of the initial and final 
state of the nucleus, angular momentum taken away by the 


183 
184 
185 
186 
187 
188 
189 
190 


19 


192 


193 


194 


195 


196 


197 


198 


199 


200 


20 


202 
203 
204 
205 
206 
207 
208 


209 


210 


21 


212 
213 
214 
215 
216 
217 
218 
219 


220 


emitted two protons, respectively. The fifth column is the 
experimental half-lives of 2p radioactivity of excited state. 
From the sixth to ninth columns, they represent the logarith- 
mic form of 2p half-lives obtained by CPPM, ELDM, GLDM 
and UFM, respectively. In general, from this table, it is ob- 
viously that the theoretical half-lives obtained by CPPM are 
highly consistent with other theoretical models, except for the 
nuclei ??Mg* (with the Q2, = 6.11 MeV) and 7°S* (with the 
Q2p = 5.12 MeV). In order to intuitively explain this phe- 
nomenon, we plot the interaction potential curve of 77Mg* 
and 79S* in Fig. 1. In Sec. IIA, the half-lives of 2p radioac- 
tivity are depended on penetration probability P which is ob- 
tained by Eq. 3 and conditions V (rin) = V (rout) = Q2p- 

Unfortunately, in Fig. 1, the whole interaction potential 
curve of ??Mg* and ?°S* can’t satisfied the conditions men- 
tioned above, so that we can’t calculate corresponding P, nat- 
urally, the half-lives can’t be obtained. 

For the reason of this phenomenon, we turn our attention 
to remarkable deformation effect. In fact, the half-lives of 
2p radioactivity are highly sensitive to proton-proton interac- 
tion due to the paring effect of valence protons. The quasi- 
classical ?He model can’t account for the experimentally ob- 
served proton-proton correlations, which indicate back-to- 
back proton emission. Moreover, for the nucleus ?7Mg* and 
29S*, their spin-parity of the initial state of the parent nuclei 
are uncertain, and the the values of Q2p are floating, all of 
these factors have nonnegligible impact on V (r) curve. 
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Fig. 2. (color online) Comparing the half-lives of the nuclei which 
obtained by CPPM in Table I with other theoretical models. 


In order to describe the agreement of the half-lives for the 
2p radioactivity of excited state which calculated by CPPM, 
partial theoretical results are plotted in Fig. 2. In this fig- 
ure, the logarithmic half-lives of the 2p radioactivity of ex- 
cited state which obtained by ELDM, GLDM and UFM are 
compared with our calculations [84, 85]. For the 140*, in 
Table 1, its experimental data are given as logioT\/2 > - 
16.12 s, the calculated half-lives by CPPM are very close 
to those from ELDM, GLDM and UFM, simultaneously the 
calculated results are better than from the some empirical for- 
mulas such as R-matrix theory [42], Sreeja formula and Liu 
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formula [85], whose half-lives are -18.12 s, -19.94 s and - 
16.85 s, respectively. In order to intuitively discuss the in- 
fluence of angular momentum £ and release energy Q2p for 
the half-lives of the 2p radioactivity from excited state, we 
select '4O* and °*Ag* as an example to analyze their contri- 
bution to the half-lives. In Table 1, it is obviously find that 
the half-lives of '4O* obtained by CPPM for the same value 
of £ nearly have three magnitudes deviation with different re- 
leased energy Q2p value (Q2p= 1.20 MeV and Q2p= 3.15 
MeV). Simultaneously, for the Yt Ag, identical Qə2p but dif- 
ferent £ also have three magnitudes while the value of angular 
momentum are 6 to 10. In fact, this conclusion is same with 
our previous work [87]. It is obviously found that either Q2p 
or l has a nonnegligible contribution to the '*O* and °*Ag* 
for their corresponding theoretical half-lives within CPPM, 
ELDM, and GLDM. The reason for this phenomenon is the 
theory of these models are similar, the penetration probatbil- 
ity P are both calculated by WKB approximation. In other 
word, these theoretical models have a strong dependency with 
Q2p and £, and we suspect the corresponding experimental 
data are not precise enough. Configuration mixing and three- 
body dynamics are commonly utilized in nuclear and molec- 
ular physics to describe systems with numerous interacting 
particles. However, due to the limitations of the phenomeno- 
logical theoretical model, which fails to account for configu- 
ration mixing and three-body dynamics, our calculations treat 
the two emitted protons as strongly correlated, thus reducing 
them to a two-body problem. As a result, CPPM cannot pro- 
vide comprehensive information on particle dynamics. Up to 
now, the experimental 2p radioactivity data of excited state 
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caused by the angular momenta are not available because the 
spin-parity of the initial state of the parent nuclei has not de- 
termined. In addition, the uncertain values of @2, also pro- 
vided strong influence on calculated results. Simultaneously, 
it is found that the half-lives of excited 2p radioactivity have 
a strong relationship with @»2, and £. Comparing with the 
theoretical results obtained by ELDM, GLDM and UFM, the 
half-lives of the excited 2p radioactivity by CPPM are reli- 
able, it maybe as a positive guideline for future experiments. 
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